Many-body Green's function study of coumarins for dye-sensitized solar cells 
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We study within the many-body Green's function GW and Bethe-Salpeter formalisms the excita- 
tion energies of several coumarin dyes proposed as an efficient alternative to ruthenium complexes 
for dye-sensitized solar cells. Due to their internal donor- accept or structure, these chromophores 
present low-lying excitations showing a strong intramolecular charge-transfer character. We show 
that combining GW and Bethe-Salpeter calculations leads to charge-transfer excitation energies 
and oscillator strengths in excellent agreement with reference range-separated functional studies or 
coupled-cluster calculations. The present results confirm the ability of this family of approaches 
to describe accurately Frenkel and charge-transfer photo-excitations in both extended and finite 
size systems without any system-dependent adjustable parameter, paving the way to the study of 
dye-sensitized semiconducting surfaces. 

PACS numbers: 71.15.Qe,78.67.-n,78.40.Me,72.40.+w 
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I. INTRODUCTION 

Promising to become a low-cost alternative to stan- 
dard silicon-based photovoltaics, dye-sensitized solar cells 
(DSSC) have been intensively studied over the past 20 
yearsJ^ The most prominent modern DSSCs, the so- 
called Gratzel cells, consist of porous layers of titanium 
dioxide nanoparticles covered by molecular dyes that ab- 
sorb sun light. While most efficient sensitizers are com- 
posed of ruthenium dye complexes, intense research is 
conducted so as to hnd molecular alternatives which are 
cheaper, easier to synthesize, and free from the resource 
limitations related to the noble metal ruthenium. As a 
promising direction, Hara and coworkers^ demonstrated 
that coumarin-based dyes, such as the so-called NKX- 
2xxx family (see Fig. [TJ, could lead to conversion effi- 
ciencies approaching that of ruthenium-based DSSCs. 

Starting from the originally tested C343 coumarin 4 
(Fig.[TJi), the introduction of (-C=C-) methine fragments 
between the coumarin unit and the terminal (-COOH) 
carboxyl group (Fig. [5J induces a red shift of the ab- 
sorption spectrum with improved light harvesting in the 
visible ranged The resulting molecular structures are rep- 
resented in Figs, l(b-d) showing the so-called (cis) con- 
formations, with (trans) structural isomers represented 
in Fig. 2 for one of them. Further, inclusion of the cyano 
(-C=N) group enhances the acceptor character of the 
combined (-COOH) and (-C=N) cyanoacrylic acid group, 
increasing the charge-transfer (CT) character of the in- 
ternal excitations, a factor believed to favor the injection 
of the photoelectron to the inorganic structure through 
the anchoring (-COOH) carboxylic unit4 Finally, the re- 
placement of the methine spacer by thiophenc chains 5 
reduces the adverse aggregation of dyes onto the Ti02 
surface, leading to the NKX-2677 dye (Fig. le) with a 
solar-energy-to-electricity conversion efficiency of 7.7%. 

Due to its donor-acceptor structure, this family of 
molecules became also a benchmark for theoretical stud- 
ies aiming at solving the well-known problem of describ- 



ing CT excitations within time-dependent density func- 
tional theory (TDDFT)^ Such difficulties arise from 
the lack of long-range electron-hole interaction when 
(semi)local kernels are used, as a signature of the miss- 
ing non-local (screened) exchange contribution^^ In par- 
ticular, TDDFT calculations with local^ hybridii^ 
and range-separated hybrid (RSH)ii~— kernels were con- 
ducted and compared to reference quantum chemistry 
coupled-cluster CC2 calculations^ in order to assess the 
accuracy of the various approaches. The central role 
of CT excitations in organic and hybrid solar cells^ 
certainly explains such efforts to develop computational 
techniques able to describe both Frenkel and CT exci- 
tations in a large variety of finite and extended systems 
with various screening environments. 

In a few recent studies, many-body Green's func- 
tion techniques, namely the so-called GW—~— and 
Bethe-Salpeter (BSE)! 9 .^^^ formalisms, have been 
used to investigate the optical excitations in a family 
of small donor/acceptor organic dyads^Zr— for which 
gas phase experiments^ are available. The low-lying 
CT excitations were found to come within 0.1 eV 
as compared to experiment] 28 ' 29 Besides these systems 
for which experimental data exist, optical excitations 
in a paradigmatic fullercnc/thiophcnc derivative accep- 
tor/donor complex^ and intramolecular CT excitations 
in a model dipeptide^ have also been studied within 
the BSE formalism. This family of Green's function 
many-body techniques has been applied since the mid- 
cightie o 17 ' 18 ' 21 at the ab initio level to study extended 
semiconductors and metals, but its applicability to the 
crucial problem of CT excitations in organic systems re- 
mains in its infancy and extensive tests and studies are 
still needed to rationalize its merits and limitations. 

In the present work, we analyze the (singlet) excitation 
energies of various coumarin-based molecules such as the 
parent C343 dye and the related NKX-2388, NKX-2311, 
NKX-2586, and NKX-2677 structures for which refer- 
ence calculations are available. We demonstrate that the 
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a) C343 




FIG. 1: (Color online) Symbolic representation of the studied 
coumarins: (a) parent C343, (b) NKX-2388 (cis), (c) NKX- 
2311 (cis), NKX-2586 (cis) and NKX-2677. The difference 
with the corresponding (trans) structures is represented in 
Fig. 2 for the NKX-2311 case. Black, white, red, blue and 
yellow atoms represent carbon, hydrogen, oxygen, nitrogen 
and sulfur, respectively. 
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FIG. 2: (Color online) Symbolic representation of the NKX- 
2311 (trans) coumarin dye. This structure can be compared 
to the NKX-2311 (cis) structure in Fig. 1(c) which differs from 
the orientation of the tv conjugated polymethine (-HC=CH-) 
bridge (orange frame), the cyano (-C=N) (blue dot-dashed 
box) and the carboxylic (-COOH) (full red box) groups. The 
donor coumarin parent is on the left side in the green dot- 
dashed frame. The anchoring group to the semiconducting 
Ti02 surface is the carboxylic unit. 



GW /BSE formalism yields the low- lying charge-transfer 
excitation energies in excellent agreement with opti- 
mized (parametrized) range-separated hybrid function- 
als and quantum chemistry CC2 coupled-cluster calcula- 
tions with a mean absolute error smaller than 0.06 eV. 
Similarly, the associated oscillator strengths are found to 
agree very well with ajusted long-range corrected func- 
tional calculations. These results consolidate the few re- 
cent evidences showing that GW /BSE calculations can 
describe charge-transfer excitations in gas phase organic 
systems with an accuracy of the order of a tenth of an 
eV without any adjustable parameter. 



II. TECHNICAL DETAILS 

Our calculations are performed with the Fiesta 
packag o 28 ' 32 " — which exploits atom-centered auxiliary 
Gaussian-basis for developing the needed two-body oper- 
ators such as the dynamical (w-frequency dependent) sus- 
ceptibilities, the screened Coulomb potential W(r, r'; cu), 
and the self-energy in the GW approximation: 

£ Gly (r,r» = f rf W ' e ,u ' 0+ G(r,r>W)^(r,r'; W '), 

277 J 

where G is the time-ordered one-particle Green's func- 
tion. The self-energy is non-local in space and time 
(energy-dependent) and accounts for exchange and cor- 
relation in the present formalism. Our auxiliary basis 
contains six exp{— ar 2 ) Gaussians for the radial part of 
each (s, p, d)-channel, with an even tempered 3 ^ distribu- 
tion of the localization coefficients a ranging from 0.1 to 
3.2 a.u. except for hydrogen where the range is set to 0.1 
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a.u. to 1.5 a.u. As such, our auxiliary basis contains typi- 
cally 54 orbitals per atom. Such a basis derives from pre- 
vious studie s 28 ' 32 ' 33 but with additional diffuse orbitals. 
Convergency tests can be found in Ref. with more de- 
tails about the present GW implementation which is sim- 
ilar to that of Rohlfing and coworkers i 36 ' 37 We however 
go beyond the plasmon-pole approximation for obtaining 
the dynamical correlations by using contour deformation 
techniques 

The single-particle states needed to build the start- 
ing screened-Coulomb potential W and Green's func- 
tion G are the Kohn-Sham DFT eigenstates in the 
local density approximation (LDA) as provided by 
the Siesta package 3 ^ with a large triple-zeta plus 
double polarization basis (TZDP) for the valence 
orbitals4^ Using such running parameters, comparison 
with GW/BSE calculations performed with a reference 
planewave code^i showed excellent agreement in the case 
of acene/tetracyanoethylene (TCNE) donor /acceptor 
complexes^ demonstrating the accuracy of the present 
Gaussian-based formalism. For sake of illustration, 
the Kohn-Sham and auxiliary basis associated with the 
NKX-2677 molecule (67 atoms) contain 1110 and 3618 
orbitals, respectively. Our structures are relaxed at the 
all-electron B3LYP 6-311G(d,p) level.i 2 - 

We start with a GW calculation designed to ob- 
tain the correct quasiparticle spectrum, namely the 
correct occupied and unoccupied single-particle energy 
levels. As observed in several studies on gas phase 
(isolated) molecular systems , 32 ' 33 ' 43 ' 45 " — the standard 
"single-shot" Go Wo (LDA) scheme tends to produce too 
small energy gaps between the highest occupied (HOMO) 
and the lowest unoccupied (LUMO) molecular orbitals. 
In the case of CT excitations with a large weight on 
the HOMO-LUMO transition, this has been shown to 
result in Go Wo/BSE excitation energies too small as 
compared to experiment, with a mean absolute error 
of about 0.7 eV in the case of the small TCNE/acenes 
complexes^ Such a problem has been largely cured 
by several group s 28 ' 29 ' 32 ' 45 " — thanks to a partial self- 
consistent cycle where the quasiparticle energies are rein- 
jected to build updated time-ordered Green's function G 
and screened Coulomb potential W, while keeping the 
eigenstates (Kohn-Sham one particle orbitals) frozen^ 4 - 
This simple self-consistent loop allows to obtain much 
better quasiparticl o 32 ' 33 ' 46 ' 47 and excitatio n 28 ' 29 energies 
as compared to experiment. Further, the observed sen- 
sitivity of the Go Wo results to the choice of the start- 
ing zeroth-order eigenstate s 28 ' 32 ' 44 ' 46 " — is significantly 
reduced^ It is such an approach that we adopt in the 
present study. 

In a second step, the neutral (optical) excitation 
spectrum is constructed by using the Bcthc-Salpctcr 
equations (BSE)^ 2 ."— that account for the electron-hole 
(excitonic) interactions. Deriving from Green's func- 
tion many-body perturbation theory, the Bcthc-Salpctcr 
Hamiltonian H BSE is most commonly written in the 
two-body product basis 0j(r)0j(r') of the Kohn-Sham 



orbitals r 4 - resulting in the following generalized eigen- 
value problem: 
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The first diagonal block: A BSE = {H dla 9 + H direct + 
jjexch^ accounts for resonant transitions from occupied 
(indexes v,v') to unoccupied (indexes c,c') eigenstates, 
with the following contributions: 
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The second diagonal block (—A*) accounts for non- 
resonant transitions from unoccupied to occupied lev- 
els. The off-diagonal C term couples resonant and non- 
resonant transitions, with C also composed of a direct 
and an exchange term which can be obtained from their 
A-block analogs by switching the (v',c') indexes. While 
all empty states are included in the construction of the 
independent-electron susceptibilities and Green's func- 
tion needed for the GW calculations, the number of 
empty states included in the BSE Hamiltonian has been 
set to 160. Convergency tests for the largest NKX-2677 
molecule indicate that increasing this number to 200 de- 
creases the excitation energy by less than 10 meV. 

In what follows, we go beyond the Tamm-Dancoff 
approximation (TDA) by mixing resonant and anti- 
resonant contributions. Namely, we do calculate the 
coupling (C) and (— C*) blocks. As shown in recent 
studies; 37 ' 55 the TDA tends to overestimate by a few 
tenths of an eV the excitations energies in small size sys- 
tems where single-particle and collective excitations can 
strongly mix. This is what we observe here with the low- 
lying excitation energies which are blue-shifted by up to 
0.25 eV for the C343 dye, and 0.18 eV in the NKX-2677 
case, when the off-diagonal coupling blocks are neglected. 



III. RESULTS AND ANALYSIS 

Our results are compiled in Table U and in Fig. [3] with 
a comparison to previously published TDDFT and CC2 
calculations. As expected, the TD-LDA and TD-PBE 
values from Ref. [l^ evidence an underestimation of the 
transition energies related to the CT character of the ex- 
citations. The mean absolute error (MAE) averaged over 
the available data points amounts to 0.47 eV and 0.44 eV 
as compared to CC2 within LDA and PBE, respectively. 
Further, the oscillator strength seems to be significantly 
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. — ttT\ 

LDA/PBE (a) 


B3LYP (i,) 


LC-BLYP (C) 


7T7V 

BNL (d) 


J1/J2 


CC2 


w 


GW-BSE 


C343 


2.96/3.0 (0.36/) 


3.32 


(0.60) 


3.36 


(0.57) 


3.5/3.4 


(0.7/0.6) 


3.44 


(0.74) 


3.44 


(0.57) 


NKX-2388 s-trans 




2.90 


(0.94) 


3.01 


(0.88) 


3.1/2.9 


(1.0/0.9) 


2.99 


(1.06) 


3.04 


(0.88) 


NKX-2388 s-cis 




2.78 


(0.87) 


2.85 


(0.80) 


2.9/2.8 


(0.9/0.9) 


2.80 


(1.00) 


2.85 


(0.80) 


NKX-2311 s-trans 




2.70 


(1.35) 


2.91 


(1.34) 


2.9/2.8 


(1.6/1.5) 


2.89 


(1.51) 


2.88 


(1.37) 


NKX-2311 s-cis 


2.35/2.35 (1.05/) 


2.56 


(1.19) 


2.73 


(1.12) 


2.7/2.6 


(1.3/1.2) 


2.71 


(1.33) 


2.67 


(1.13) 


NKX-2586 s-trans 




2.50 


(1.71) 


2.81 


(1.83) 


2.8/2.6 


(2.1/2.0) 


2.81 


(2.01) 


2.74 


(1.88) 


NKX-2586 s-cis 


2.10/2.15 (1.23/) 


2.40 


(1.55) 


2.66 


(1.52) 


2.6/2.5 


(1.7/1.7) 


2.66 


(1.74) 


2.57 


(1.59) 


NKX-2677 




2.23 


(1.49) 


2.67 


(1.76) 


2.7/2.5 


(2.0/1.8) 


2.71 


(2.17) 


2.56 


(1.69) 



TABLE I: Calculated lowest (So Si) singlet transition energies (in eV). The GW-BSE results calculated in the pre sent 
study are compared to the TD-LDA, TD-PBE, TD-B3LYP, TD-LC-BLYP, TD-BNL and CC2 calculations from Refs. tuj4jj. 
We provide the B3LYP results from Ref. [Tl| which are in excellent agreement with the results of Ref. [Isl . The numbers in 
parenthesis indicate the associated oscillator strengths. 
"Ref.llC. 
b Ref. U 
c Ref. 12. 
d Ref. " . 



underestimated, in particular, surprisingly, in the case of 
the smaller C343 parent molecule expected to show the 
smallest CT character. 

Thanks to its 20% of exact exchange, the B3LYP re- 
sults of Ref. [ll] (filled black circles in Fig. [3]) indicate a 
reduced MAE of 0.2 eV. However, the discrepancy can 
still be as large as 0.48 eV for the NKX-2677 structure. 
This is certainly the signature that in the long-range 
charge-separation limit, the restricted amount of exact 
exchange in the B3LYP functional 5 ^ is not enough to ac- 
count for the correct electron-hole interaction. To illus- 
trate that point, we plot in Fig. [4] the Kohn-Sham (a,d) 
HOMO and (b,e) LUMO cigenstates associated with the 
C343 and NKX-2677 structures. While the HOMO states 
arc found to be rather delocalized, the LUMO in the 
NKX-2677 dye is clearly much more localized close to 
the electron-acceptor cyanoacrylic group, resulting in an 
enhanced CT character as compared to the C343 parent 
molecule. 

The nature of the transitions can be better quantified 
by studying the electron and hole spatial localization in 
the excited states. This can be achieved by taking the 
expectation value of the electron/hole position operator 
S(r — r e /vh) over the two-body ^(r e ,r),.) BSE excitonic 
wavefunction, leading to an electron/hole probability of 
presence averaged over the hole/electron position. The 
resulting densities are provided in Fig. |U[c,f) with an 
isocontour representation for the C343 and NKX-2677 
molecules. These densities allow to obtain the mean elec- 
tron/hole positions (see red arrows in Fig. @| and the 
related average electron-hole separation distance which 
amounts to 3.2 A in C343. This clear CT character is 
certainly at the origin of the difficulties met by LDA or 
PBE to describe such an excitation. In the NKX-2677 
case, this average distance increases to 4.6 A as a signa- 
ture of the enhanced CT character, explaining that the 
B3LYP results significantly worsen from C343 to NKX- 
2677. 

We now come to the central results of the present 
study, namely the many-body perturbation theory data. 



In contrast with the LDA, PBE or even B3LYP results, 
our GW-BSE values (empty red squares in Fig. are in 
much better agreement with the CC2 data points, with 
a mean absolute error of 0.06 eV. Such an agreement 
is remarkable accounting for the fact that the present 
GW /BSE approach does not contain any adjustable pa- 
rameter. Concerning the longest NKX-2677 dye, which 
shows the largest discrepancy with CC2 calculations, we 
observe that our result falls within the values provided 
by the range-separated hybrid 5 -^ (RSH) BNL functional 
study^ where two different strategics to optimize ab 
initio (non-empirically) the range-separation parameter 
have been tested^ As compared to the RSH-BNL study, 
our GM^-BSE results differ by a MAE ranging from 0.04 
eV to 0.07 eV, that is well within 0.1 eV. As empha- 
sized in Ref. [lj, the CC2 approach is also not free from 
approximations 59 and differences of the order of 0.1 eV 
as compared to more accurate e.g. CASPT2 calculations 
are certainly to be expected. 

Clearly, as compiled in Tabic HI TDDFT calculations 
with the LC-BLYP functional^ also provide excellent 
results^ with a 0.03 eV MAE as compared to CC2, 
smaller than our GW^-BSE 0.06 eV MAE value. However, 
as emphasized in Ref. [l2l . the range-separation parameter 
[i in the LC-BLYP study has been precisely adjusted to 
minimize the root mean square error with CC2 calcula- 
tions. The best-fit /i value for these systems (/x=0.17 a.u.) 
is found to be much smaller than the original /x=0.33 ad- 
vocated by Iikura and coorkers^. The strong dependence 
of the excitation energies as a function of fj. indicates that 
the choice of the originally recommended /i=0.33 value 
would lead to a significant overestimation of the tran- 
sition energies (by as much as 0.3-0.4 eV, see Fig. 4 of 
Ref. [l2l) . This leads to the standard question of the choice 
and transferability of the range-separation paramcter(s). 
One observes however that with the best-fit [i value, the 
correlation between LC-BLYP and CC2 results is very 
remarkable, showing that this class of systems can be 
described by a unique parameter. 

Bearing important consequences on the use of such 
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2.6 2.8 3 3.2 3.4 

CC2 Excitation energies (eV) 

FIG. 3: (Color online) Calculated lowest singlet excitation 
energies (in eV) as a function of the coupled-cluster CC2 ref- 
erence values. Results in perfect agreement with the CC2 
calculations should fall on the first diagonal (black line) . The 
present GW-BSE calculations (empty red squares) are com- 
pared to the TD-B3LYP results of Refs.[H|GJ (filled black cir- 
cles) and the "adjusted" (see text) TD-LC-BLYP data from 
Ref. Q3 (blue triangles up). The coumarins name are indi- 
cated by their number (removing the NKX prefix) with (c) 
standing for -cis and (t) for -trans. The axes physical length 
is scaled according to their respective energy range. 



dyes in DSSCs, our GW-BSE calculations confirm the 
range-separated hybrid TDDFT and CC2 data leading 
to the conclusion that the onset of absorption is signifi- 
cantly red shifted with increasing size length. This evo- 
lution is in clear contrast with the behavior of CT exci- 
tations in well separated gas phase donor/acceptor dyads 
where the exciton binding energy scales as the inverse dis- 
tance between the two molecules, leading to an increase 
of the absorption energy onset. However, contrary to 
well separated donor /acceptor systems, the quasiparticle 
HOMO-LUMO gap in donor /acceptor dyads connected 
by a conducting 7r-conjugated bridge does not remain 
constant with varying bridge length. This is clearly ex- 
emplified in Table IJ where the GW HOMO-LUMO gap 
is found to quickly decrease from the C343 molecule to 



C343 Nkx2677 




FIG. 4: (Color online) (a) and (b) Isocontour representation 
of the C343 HOMO and LUMO Kohn-Sham eigenstates. Dif- 
ferent colors indicate different signs of the wavefunction. In 
(d) and (e), similar plots for NKX-2677. In (c) and (f), iso- 
contour representation of the electron (yellow) and hole (light 
blue) probability distribution for the lowest C343 and NKX- 
2677 singlet excited states, respectively, as obtained within 
BSE. The red arrows indicate the average hole (left arrow) 
and electron (right arrow) positions. The C343 and NKX- 
2677 molecules are not represented on the same scale. 



the longest NKX-2677 dye. Except for the large variation 
of the electronic affinity (EA) from the C343 parent to 
the NKX-2388 system, analysis of the GW HOMO and 
LUMO quasiparticle energies indicates that this gap re- 
duction stems both from a destabilisation of the HOMO 
(decrease of the ionization potential IP) and a stabiliza- 
tion of the LUMO (increase of the EA). Subtracting the 
GW HOMO-LUMO gap from the singlet GW-BSE ex- 
citation energy, one finds that the electron-hole (exci- 
tonic) binding energy Eb decreases with donor-acceptor 
distance, but this effect is not strong enough to coun- 
terbalance the decrease of the HOMO-LUMO gap. This 
is an important feature which explains that the longest 
NKX-2677 dye is most efficient in harvesting photons in 
the visible range. 



We close this study by commenting on the oscillator 
strength associated with the calculated transitions (see 
the numbers in parenthesis in Table U and Fig. [5]). Our 
GJU-BSE values (red empty squares) lead to an excel- 
lent agreement with the parametrized LC-BLYP (blue 
up triangles) calculations, with a perfect match for the 
C343 and NKX-2388 molecules. Such an agreement is 
remarkable given the fact that the two approaches are 
very different. Both GIU-BSE and LC-BLYP provide 
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IP 


EA 


Gap 


E B 


C343 


7.21 


0.60 


6.61 


3.17 


NKX-2388 s-trans 


7.12 


1.33 


5.79 


2.75 


NKX-2388 s-cis 


7.11 


1.50 


5.61 


2.76 


NKX-2311 s-trans 


6.93 


1.55 


5.38 


2.50 


NKX-2311 s-cis 


6.92 


1.71 


5.21 


2.54 


NKX-2586 s-trans 


6.76 


1.70 


5.06 


2.32 


NKX-2586 s-cis 


6.77 


1.85 


4.91 


2.34 


NKX-2677 


6.52 


1.62 


4.90 


2.33 



ysis is needed to understand these variations from one 
type of calculation to another. Beyond differences, one 
should emphasize that all studies agree on the fact that 
the NKX-2677 dye presents one of the largest oscillator 
strength, yet another factor explaining that it leads to 
one of the largest conversion efficiency in the coumarin 
family. 



TABLE II: Calculated GW ionization potential (IP), elec- 
tronic affinity (EA), and HOMO-LUMO gap. The exciton 
binding energy Eb is the GW HOMO-LUMO gap minus the 
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FIG. 5: (Color online) Theoretical strength as a function 
of the GW-BSE excitation energies (in eV). The BSE val- 
ues (red empty squares) are compared to the CC2 (black di- 
amond), LC-BLYP (blue triangle up), and the BNL J1/J2 
(empty green circles) results. Results for the NKX-2677 and 
NKX-2586 (s-cis) in the dashed box are reproduced in the 
up-right inset. 



lower values as compared to CC2 calculations, the BNL 
data yielding somehow intermediate results. The largest 
discrepancies occur for the NKX-2677 molecule, with a 
22% difference between the GW-BSE and CC2 values. 
For this system, GW-BSF, and LC-BLYP agree again ex- 
tremely well, within 4%, while the two BNL values (see 
Inset) show their largest spread, indicating that the os- 
cillator strength for this structure is very sensitive to the 
chosen formalism and related parameters. Further anal- 



IV. CONCLUSION 

In conclusion, we have studied within the many-body 
Green's function GW and BSE perturbation theory the 
excitation energies of a family of coumarin dyes recently 
shown to be very promising candidates for replacing 
ruthenium-based chromophores in dye-sensitized solar 
cells (DSSC). In such donor-bridge-acceptor molecules, 
the lowest singlet excitations are characterized by a 
charge-transfer character that varies with the length of 
the 7r-conjugated bridge. As a result, TD-B3LYP calcu- 
lations can lead to an error as large as ~0.5 eV as com- 
pared to reference quantum chemistry coupled-cluster 
CC2 calculations, despite the 20% of exact exchange con- 
tained in its functional form. We demonstrate that the 
GIE/BSE approach leads to an excellent agreement with 
CC2 data with a mean absolute error of the order of 
0.06 eV for the excitation energies. Such an accuracy 
is comparable to the best results provided by TDDFT 
calculations with optimized long-range corrected range- 
separated hybrids, but with a parameter-free approach 
that performs equally well for extended insulating or 
metallic systems and gas phase organic molecules. Such 
an excellent agreement is also demonstrated for the re- 
lated oscillator strengths. The ability of the GIK/BSE 
approach to describe both Frenkel and charge-transfer 
excitations in finite size molecular systems or extended 
semiconductors originates in particular from the use of 
the screened Coulomb potential W that automatically 
adjusts the strength and range of the Coulomb interac- 
tions. This flexibility may proove as a significant advan- 
tage in the study of DSSCs where both the organic dye 
and the extended TiC>2 semiconductor must be treated 
with sufficient accuracy. 
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